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ABSTRACT 


Doubly differential cross sections (ddcs) for 
ejection of secondary electrons of various energies 
at various angles from helium bombarded by 200- 
and 300-keV protons have been calculated and com- 
pared with the recent experimental results of Rudd, 
Sautter, and Bailey. The wave functions employed 
by Mapleton to calculate ionization cross sections 
are used in the Born approximation in a first calcula- 
tion. In a second calculation, the Roothaan four- 
term wave function is used for the ground state of 
helium. The two calculated results differ very little. 
General agreement between theory and experiment 
for the ddcs is found, but discrepancies are noted 
in the forward and backward directions. Changing 
the ground- state wave function from the hydrogen- 
like to the Roothaan four-term wave function has 
very little effect on the ddcs. 
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SUMMARY 


Doubly differential cross sections for ejection of secondary electrons of 
various energies at various angles from helium bombarded by 200- and 300-keV 
protons have been calculated and compared with the recent experimental 
results of Rudd, Sautter, and Bailey. The wave functions employed by 
Mapleton to calculate ionization cross sections are used in the Born approxima- 
tion in a first calculation. In a second calculation, the Roothaan four-term 
wave function is used for the ground state of helium. The two calculated 
results differ very little. General agreement between theory and experiment 
for the cross sections is found, but discrepancies are noted in the forward 
and backward directions. 


INTRODUCTION 


The cross sections for ejection of secondary electrons (with various 
energies and at various angles) from helium bombarded by 50- , 100- , and 
150-keV protons have been determined theoretically and reported (ref. 1). 
From these cross sections, which may be termed doubly differential cross 
sections (ddcs), cross sections differential in ejection energy may be obtained 
by integrating the ddcs over the ejection angle. Cross sections for ionization 
are then obtained by integrating the differential cross sections over the energy 
of ejection. Several authors - Erskine (ref. 2), Dalgarno and McDowell 
(ref. 3), and Mapleton (ref. 4) - have published calculated values for the 
ionization cross sections, but, before the results given in reference 1 
appeared, no theoretical results were available for the ddcs for comparison 
with experimentally determined ddcs. The Born approximation was utilized 
in reference 1, and products of hydrogen-like wave functions were used to 
approximate the wave functions of helium in both the ground state and the 



continuum states. The results of these calculations were compared with the 
experimental results of Rudd and Jorgensen (ref. 5). At the time reference 1 
was published, the experimental results were given only at incident proton 
energies of 50-, 100-, and 150-keV. In these cases, the ddcs agreed in a 
general way with the experimental results, but discrepancies were noted, 
particularly in the forward and backward directions. Now experimental data 
at incident proton energies of 200- and 300-keV are available (ref. 6), and the 
Born approximation is more accurate at these levels. 

The two purposes of this paper are : 

1. To perform the calculations shown in reference 1 at the higher 
incident proton energies and to compare the results with the experimentally 
determined values. 

2. To give the results of calculation of the ddcs using the four-term 
Roothaan ground- state wave function (ref. 7) for helium instead of the 
hydrogen-like wave function used in reference 1. The Roothaan wave function 
is determined by self-consistent field methods and agrees very well with the 
Hartree-Fock wave function. 

The motivation for purpose (1) is that, since the Born approximation is 
a high-energy approximation, the results from (1) should agree with the 
experiment better than did the results given in reference 1. The reason for 
(2) is to determine the sensitivity of the ddcs to the choice of a ground- state 
wave function. 


SYMBOLS 


A momentum change vector of incident protons 

a^ radius of first Bohr orbit 




defined by equation (19) 



E = 7r/2 - 


E. 


E 


n 



e 



energy of ejected electron 
final binding energy of helium atom 
initial binding energy of helium atom 
electronic charge 


f = 



F 


f n,k 
*0 


G 

g= 2Z 1 
X D 

J o 

A 

K 

n 

y\ 

K 0 


k 


- Z 2 ^j(z i - Z - A 2 (z i + Z 3 ) + 2ZgAk cos 6 

hypergeometric function 

Born matrix element 

defined by equation (21) 

^k^Zj - + ZgA cos ej 

Planck's constant divided by 2 tt 
doubly differential cross section 

Bessel function of order zero 

final proton wave vector 

initial proton wave vector 

vector of electron 
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k 2 + A 2 + Z^ 2 - 2Ak cos 6 
[(A 2 - k 2 + Z*f + 4k* Zl 2 ] 1/2 


a 

orbital angular momentum quantum number 

A 

n o 

incident proton direction 

A 

n i 

direction of proton scattering 

A 

n 2 

direction of electron ejection 

n' 

direction of change of momentum of incident proton 

n, s, p 

electron states 

(cos e) 

Legendre polynomials 

P 

probability of ejecting an electron into the solid angle da> with 
momentum between k and k + dk 

1 

CM 

^ ^ 
11 

« 

k 2 + Zl 2 ) 2 + 4k 2 Zl 2 

/V 

r 

radius vector 

r r r 2 

radius vectors of electrons 1 and 2 

S = k 2 + A : 

2 2 

+ Zj - 2Ak cos 0 

W 1 

defined by equation (9) 

W 2 

defined by equation (10) 

X, Z 

spherical coordinates 

Z 

nuclear charge 
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constant 


P 


r 

6,p 

77 = r(l 

e 

e k 

e 

r 

V- 


Gamma function 

spherical coordinate angles of scattered proton 
cos 6) 

angle between r and the direction of ejection 

angle between direction of ejection and momentum change vector 

angle between radius vector of ejected electron and momentum 
change vector 

reduced mass of the two-body system 


v = 



(k + A^ + Zj 2 

(k - A ) 2 + Zj 2 



Roothaan ground- state wave function of helium 
continuum wave function of hydrogen-like atom 

4>^ where i = 0 

ground- state wave function of hydrogen-like atom of nuclear 
charge Z^ 


X 


angle between momentum of ejected electron and incident proton 
direction 

positive energy function for a Coulomb field of charge -Z^e 
wave function of the continuum state 



iJ/q ground-state wave function of helium 

dfl solid angle of scattered proton 

do> solid angle of electron ejection 


THEORY 


The probability p of ejecting an initially bound electron into the solid 
angle do> with momentum between k and k + dk and scattering the incident 
proton into the solid angle dfl is 


p = f^ 1 ’ ^ 2 k 2 dk dco dfl 


( 1 ) 


n k 

where f^ ’ is the Born matrix element for the transition from the ground 

state of helium to the continuum state with the bound electron in state n and 
the free electron with momentum k. If r and r ^ are the position vectors 

of the two electrons, then, according to Mapleton (ref. 4), 


•a' ^ 

lA 


f ° ' /^ k > n ( r !’ r2 ) e ^o( r l’ r 2) dr l dr 2 ® 


where 


p = reduced mass of the two-body system 
x(/q = ground- state wave function of helium 
* 

\f/, = complex conjugate of the wave function of the continuum state 

Kj n 

/\ A /V 

A = Kq - K n , where Kq and are, respectively, the initial 

and final wave vectors of the proton in the center-of-mass 
system 
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! 


e = electronic charge = 1 

•fi = Planck' s constant divided by 2 tt 

Only final states in which the remaining bound electron is a Is electron will 
be considered; hence, n will be dropped from the notation for the matrix 
n k 

element f^ ’ . The justification for restricting the final- state bound electron 

to a Is electron is found in reference 4, where calculated results are given 
for ionization cross sections in which the final-state bound electron is in the 
Is, 2s, or 2p state. The cross sections for the latter two were found to be 
at least a factor of 100 less than the first one. Hence, ignoring the contribu- 
tion such states make to the ddcs is probably justifiable. 

For the ddcs 1^, equation (1) must be integrated over dS2, the solid 
angle of the scattered proton. 


I 


D 





(3) 


The calculations in reference 1 use the wave functions as stated in case III of 
reference 4, where the helium wave functions are approximated by products 
of normalized hydrogen-like wave functions. These wave functions are 
repeated here only for completeness. Thus, 

V) 2= * o= M Zi l r ^°( Zi l r 2) (4) 

where = 1.6875 and where is the hydrogen ground-state wave func- 
tion. 

Then 

^k,n( f l’ f 2) = WtB) j^ k (z 3 |?i) * n ( z 2 | r 2 ) + '''k^sfc) *n( Z 2 | r l)] (5 > 


I 
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where * s positive energy function for a Coulomb field of charge -Z^e, 

and (f> is the ground- state wave function of hydrogen with Z 9 = 2. To insure 
n A 

orthogonality of the initial and final states, Mapleton (ref. 4) let 


•#' k = 0 k (1|r) -^k O(1|r) + \ o ( z i| r ) 


( 6 ) 


where 


Z x = 1.6875 


0 , 


(z \r\ - 1 f in . . I 1 / 2 ex P( ikr ) 

\ 3| J 2n j^l - exp (-i27m)J r(l + n) 

x J Jq ^2(iku??)^^J exp(-u)u n du 
1 T in "1 1/2 

= sFfsnro j_~ ^(~-i2im) j exp(ikr) 


x r(C + 1 + n) (2kr)^ P. (cos 0) 

1=0 


x F(l + 1 + n, 2i + 2, -i2kr) 


(7) 


and 


n = Zj/ik 
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T] = r(l - cos 6) 


I 



where 6 is the angle between r and the direction of ejection. Here, 
is the St = 0 part of <p ^ . 


0 


If equations (4), (5), and (6) are substituted into equation (2) for the 
Born matrix element and the indicated integrations over dr ^ and dr 2 are 

carried out, the results can be written as follows : 


f 


Q k = C Q |w x - W 2 (l) + W 2 (l. 6875) 


( 8 ) 


where 




C 0 " .2 
ttA 


Z 2 =2 


1= f K( Z 3\ T ) elA ' T e 1 

/ *0, \ •7'^ -Z,i 


dr 


The results of performing the integrations over dr for and W 2 are as 
follows : 


W x = -4 


z 3/k 

1/2 , . 
f + lg 

- - pvn 

r- z 3 / z 3 

1 - exp ^-i27rZ 2 ykj 

^172 exp 

k bi + l \k lj + b, J 


( 9 ) 
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w = - — 

2 A 


z zA 


1 - exp^irZ^/k'j 


1/2 


where 


[(A 2 - Zl 2 - k 2 ) s^^Acos J , X 
£(k + A) 2 + Zj 2 (k - A) 2 + Zj 2 J ' 


f = (k 2 - - Z 3 )- A 2 (z x + Z g j + 2ZgAk cos 9 

l[ k ( Z l “ Z s) + Z 3 A COse ] 


g= 2Z 


E = 77/2 - tan -1 ^A 2 - k 2 + ^2kZ 1 


L = In 


k 2 + A 2 + Zj 2 - 2 Ak cos 9 


17.2 .2 „ 2\2 ,, 2 „ 2 

IA - k +Zj I + 4k Zj 


1/2 


S = k 2 + A 2 + Zj 2 - 2Ak cos 9 


R = (a 2 - k 2 + Z x 2 ^ 2 + 4k 2 Z 1 J 


( 10 ) 


y = W ln 


(k + A) 2 + Z x 2 
(k - A) 2 + Zj 2 
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n 


and 6 is the angle between A and k. Now, 1^ may be written as 

I D = C 0 2 k 2 j |w i -W 2 (l)+W 2 (1.6875)| 2 dfi (11) 

To carry out the double integration over dfl for 1^ , the reference 

direction must be changed to a fixed direction in space. In this case, using 
the incident proton direction as a reference direction is most convenient. 
Suppose the electron is ejected in the direction (X,0) and the proton is scat- 
tered in the direction (6, p), as shown in figure 1. If Eq and E n are the 

positive initial and final binding energies of the helium atom, then 

A= (2 n) 1/2 1 2E 0 - E' - 2^E 0 - E')E Q j 1/2 cos «| 1/2 (12) 

where Eq is the incident energy, and E 1 = E.. + Eq - E n where E.. is the 
ejection energy of the electron. Also 


cos 6 = 


(2p) 


1/2 


( K 0 ' K n cos 6 ) 


(2p) 


1/2 


cos 


X 


+ 




sin 6 sin X cos p 


(13) 


After using the conservation of energy and the binomial expansion in terms 
of the small quantity a = E' /Eq, A can be written 


A = ^2pEq^^ 2 ^2(2 - a) sin 2 (6/2) + cos 60 a 2 + ^ a 3 + . . .^j 


1/2 


(14) 
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and cos 6 can be written 



The double integration indicated in equation (3) can now be carried out 
over 6 and p. A computer program was written to perform these integra- 
tions for input parameters of E Q , E. , and X. Samples of results of these 

3 

calculations are presented in figures 2 through 5, which also show the experi- 
mental results for comparison. 

The calculation of the ddcs using the Roothaan ground- state wave function 
can be carried out with the procedure used above, except that r^j in 

equation (4) is replaced by the Roothaan ground- state wave function: 


-ff-r -a^r -o' r -o' r 

4»o = + fi^re + /3 g e + fi^re (16) 

with 

jSj = 0. 429299 
fin = 0. 110917 
fi 3 = 1. 47974 
fi 4 = -0. 169854 
= 3. 0 

“ 2 = 14 

12 


However, this change greatly increases the work and the numerical computa- 
tion. The Born matrix element is now 


f ° k = ^ / K ( z s| r ) elA-r v r > df < i7 > 

where 

C= /* o( r 2)*( r 2) df 2 

and 0(r) is the ground-state wave function of hydrogen with Z = 2. 

Since ^(r) is the sum of four terms, the Born matrix element in 
equation (17) can be written as the sum of four integrals. From observation 
of the form 4>Q(r), f^ can then be written in terms of of equation (9) 
and aW 1/ /9Z 1 . 


k 4uC 9W 1 0W 1 

f o = ^i* 3 ’ « - ^2 -0Z7 (3 ’ « + W 1 * 4 ’ D - ^4 -az 1 4 ’ D 

A 1 1 


(18) 


Here, W^Z^, Z^j refers to equation (9) with considered as a function 
of the parameters Z^ and Z^. In the first two terms of equation (18), 

Z i = a i , and in the second two terms Z ^ • Hence, most of the integra- 

k 

tion for f^ has already been carried out. The derivative of with 
respect to Z ^ is tedious, but straightforward, to carry out. Introducing 


13 



( 19 ) 


C = -4 
w 


V k 


i - 


exp ^- 27 TZ^y^kj 


1/2 


then gives 


0W, 


— = W 

aZj l 


>3 8L ( Z l) (. Z sV E ( Z l) 

1 k aZj + V " k/ az x 


- Z 3 E ( z l)/ k 

s 3 (Zi)r 3 / 2 ( z i) 


C exp 
w ^ 


x exp 



_~F L ( Z l) + E ( Z l)J S ( Z l) R ( Z l) 


af 


( z i) 


az. 


, t 8g ( z i) 

+ az. 


- G 


( z i)l 


( 20 ) 


with 


G ( z iH f ( z i) +ig ( z i)] 


4z i r ( z i) 


s( z i) 3R( Zl )' 


az. 


( 21 ) 


The noted derivatives are as follows : 


2 2 2 

-UZ = k 2 -A 2 -3Z 1 2 + 2Z 1 Z 3 


( 22 ) 


*( z i) jN 


3Z. 


Z. +2kZ l 


(23) 
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( 24 ) 




= 2Z 


1 


(25) 




= 4Z 





(26) 


(27) 


Now, Ip can be found by carrying out the double integration over the proton 

scattering angle as before. Some results are shown in figures 2 through 5 for 
comparison with the experiment and also with the first calculation. 


RESULTS AND DISCUSSION 


The theoretical results are in general agreement with the experimental 
results (see figs. 2 through 5), but discrepancies occur in the forward and 
backward directions. Also noted is that the agreement between theory and 
experiment improves for ejection angles less than 90° as the ejection energy 
increases, but agreement seems to worsen for angles greater than 90°. It is 
probably noteworthy that the results for 300-keV incident protons are neither 
significantly better nor worse than the same results for 150-keV protons as 
given in reference 1, except for one aspect: for a fixed ejection energy, the 
theory predicts a maximum value of the ddcs (considered as a function of ejec- 
tion angle) whose location depends on the ejection energy and always occurs 
between 0° and 90°. The experiment always shows an absolute maximum in 
the forward direction or at 0°. Only at high incident proton energies and high 
electron ejection energies does the experiment show a relative maximum to 
match the theoretically predicted maximum. Even at these higher proton and 
electron energies, the experimental results sharply increase again to an 
absolute maximum in the forward direction, but the theory predicts decreases 


15 


in this region. The reason for this discrepancy is not known. It seems 
likely that this discrepancy is due to error in wave functions or errors 
inherent in the Born approximation. 


CONCLUDING REMARKS 


As shown by these calculations, changing the ground- state wave function 
from the hydrogen-like to the Roothaan four-term wave function has very 
little effect on the ddcs. This is disappointing, considering the extra labor 
involved in using the Roothaan ground- state wave function instead of the 
hydrogen-like wave function. However, it does indicate that for future 
calculations the simpler wave function may be used. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, July 15, 1966 
030-00-00-CF-72 
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n o direction of incidence; nl direction 
of scattering; n2 direction of ejection; 
n' direction of change of momentum of 
incident proton 


Figure 1. - Directions of the incident and scattered particles 




f 




Cosine ol ejection angle 

Figure 2. - Differential cross section for production of secondary electrons by 

200- keV protons I . Figure 3. - Differential cross section for production of secondary electrons by 200-keV protons n. 




Figure 4. - Differential cross section for production of secondary electrons by 300-keV protons I, Figure 5. - Differential cross section for production of secondary electrons by 300-keV protons D. 


NASA- Langley, 19tib S-113 


19 


I 



"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington/ D.C. 20546 



